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Abstract 

I present an exactly solvable model of a pseudogap with two zero-energy 
fermion modes coupled to each other by a classical source of frequency ujq 
and strength |A|. A suitably defined fermion propagator has an infinite num- 
ber of poles at frequencies that are multiple integers of ujq. In the adiabatic 
limit, ujq <C |A|, the situation is qualitatively different from the static case 
ljq = 0: the residue of the pole at uj = (a remnant of the bare fermion) 
vanishes linearly with uiq, a result that could not be anticipated by perturba- 
tion theory; the multiple poles of the propagator coalesce into a continuum 
instead of forming two single poles at ±|A|, which should be interpreted as 
inhomogeneous broadening of the Bogoliubov quasiparticles. 
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Typeset using REVT^X 



The notion of a pseudogap is relevant to physics of fermionic systems with slowly fluc- 
tuating excitations of the Bose type. In a state with long-range order, fermions are strongly 
scattered by a condensate of these excitations, which splits the Landau quasiparticle peak 
into two and thus creates a gap in the fermion energy spectrum. In the disordered state, slow 
fluctuations can produce a qualitatively similar effect. Underdoped cuprate superconductors 
offer an example of such behavior. Without carrier doping, the cuprates are antiferromag- 
netic (AFM) insulators. Although hole doping destroys the AFM long-range order, slow 
spin fluctuations create a remnant of the Mott-Hubbard energy gap (a few tenths of an 
electronvolt [T|]). There is very little spectral weight at lower energies, which makes these 
materials rather poor conductors. 

By using a low-level approximation to estimate the fermion self-energy (a bare intermedi- 
ate fermion and a reasonably guessed particle-particle susceptibility), Kampf and Schrieffer 
have found an interesting feature in the fermion spectral function in a pseudogap regime. 



In addition to two peaks at energies ±ek = ±yek + |A| 2 (as expected in the ordered state), 

they found a small third peak at the old location ek, a remnant of the Fermi-liquid state. As 
the fluctuations slow down, the third peak gradually disappears. If true, this finding could 
give an insight into the formation of coherent spectral features (i.e., fermionic quasiparticles) 
in the cuprates upon doping. 

Deisz, Hess, and Serene have attempted to verify these results in the conserving FLEX 
approximation j|] and found no midgap peak. In our opinion, this negative result is related 
to a particular form of the self-consistently calculated spin susceptibility: while in || it 
was assumed that magnons have a typical real frequency u , the self-consistently calculated 
susceptibility clearly has a diffusing form. Whether or not propagating magnons can exist 
in the pseudogap regime is certainly an open question. Even if they can, it is far from clear 
whether the midgap band is a robust feature or an artifact of a low-level approximation. 

With this in mind, I have studied a toy model which contains the essential features of 
the Kampf-Schrieffer theory, yet permits an exact solution. Basing on the results obtained 
within this model, one can make the following conclusions. 

(a) The approximation of Kampf and Schrieffer gives reliable results in the limit of fast 
fluctuations co>o/|A| ^> 1, which corresponds to a Fermi-liquid state. 

(b) The quasiparticle peak at u = persists well into the pseudogap state (co>o/|A| C 1). 
Its spectral weight vanishes as co>o/|A|, which is a non-perturbative result. 

(c) In the adiabatic limit Uq — > 0, the two quasiparticle peaks at ±ek are strongly broad- 
ened. Broken ergodicity is necessary to restore their sharpness. 

(d) When the energies of the coupled fermion modes are not equal, the spectral function 
in the adiabatic limit is a continuum with three peaks, none of which is a pole of the fermion 
propagator. 

The toy model is defined as follows. Two fermion modes a\ and interact with each 
other via a classical, time- dependent external source whose frequency equals Uq] the Hamil- 
tonian of this system is 



By an appropriate unitary transformation, A can always be made positive, A = A* = |A|. 
It should be noted that a similar model of two fermion modes interacting with a single 
(quantum) boson mode of frequency u has been considered by Schrieffer || . 





(1) 
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If the external source is static, uiq = 0, the Hamiltonian is readily diagonalized: there 
are two eigenmodes a± = (ai ± g^/v^ with energies =f|A|. When cu ^ 0, symmetry with 
respect to time translations is broken and the time-ordered propagator matrix G a/ g(t,t') = 
—i(ip\T[a a (t)a'Jt')]\ijj) becomes a function of two time variables. In order to restore this 
symmetry (and conservation of energy with it), I perform translations t — > t + r, t' — > t' + r 
and average G(t+T, t'+r) with respect to r. This process is similar to restoring translational 
invariance in the impurity problem by averaging over impurity positions throughout the 
crystal. The resulting propagator 

G(t-t') = (G(t + T,t' + T)) T (2) 

is a function of t — t' only. Treating if as a perturbation, one can write out a series for 
G(t, t') by iterating the matrix Dyson equation 

G(t,t') = G {0) (t -t') + |A| J dt" G {0) (t - t") ai coscuot" G(t",t'). 

Upon averaging over time shifts, diagrams with an odd number of cosine factors vanish, 
making the series for G diagonal (this would not be so in the static case): 

G(t-t') = G (0) (t-t') 

+ 2 / dtldh G( ° )(t ~ t2) G( ° )(t2 ~ tl) COSUJ °( t2 ~ h ^ G(0> ^ ~ ( 3 ) 
+ . . . 

A generic term of the perturbation series for G can be readily written out, most easily for 
the Fourier transform of (0) G{u) — see Fig. |l| Note that the effective coupling constant is 
^|A| 2 /4. 

To second order in | A | , the self-energy is 

e «m = J^L + |A|2 |A|2 ^ 



It is similar to E of Kampf and Schrieffer in that it vanishes with a negative slope at uj = 0. 



Therefore, to this order, G(u) = — S(u;)] has three poles, at u — iy^o + l^l 2 /2 and 
at to = 0, the latter with the residue 

Z ° ~ u;g+ |A|V2' ( } 

which apparently vanishes as 2o;q/|A| 2 in the limit of slow fluctuations. Simultaneously, the 
other two poles are approaching ±|A|/v2 and their residues tend to 1/2, as in the static 
case. When fluctuations are fast, the dominant pole is at to = 0. Thus a picture of a smooth 
crossover from one to two peaks in the spectral function emerges. 

This smoothness is, however, misleading. As has already been noted, there are off- 
diagonal contributions to the propagator matrix in the static case, but not in the fluctuating 
case. Therefore, the adiabatic limit (ujq — > 0) is not expected to resemble the static situation 
(co>o = 0). In addition, higher powers of |A| in the self-energy can be neglected only when |A| 
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is small, which for low fermion frequencies u translates into |A| <C ujq- Another indication 
of a nonperturbative character of the adiabatic limit comes from the exact expression for 
the residue of the u = pole obtained below. As g = |A| 2 /4 — > oo, the behavior of z is 
nonanalytic: z ~ const g^ 1 ^ 2 . 

Exact solution. By making a unitary transformation, the Hamiltonian (|1|) can be 
diagonalized (in the same way as in the ojq = case): 

H = |A| cosu^ (a_a- — a) + a + ). 

Time dependence of the annihilation operators is given by 

a±(t) = a±(Q) exp ( ±z- — - smuot ) . 

\ uj J 

After choosing a state vector one can immediately write down the propagator G(t—t'). A 
particular choice of |?/>) does not affect the form of G(u) but rather determines the integration 
path in the complex plane of a;. If we select the state annihilated by a_ and a\, 



G±(t,t') = =Fzexp 
and, after averaging over time shifts, 



±i (sin uj$t — sin uot ) 

u 



6(±t =f t') 



\ wo 2 J 

where Jq(x) is a Bessel function. Note that G(t,t') and G(t — t') behave differently in the 
limit ujq —>■ 0. Without averaging, fermions a± have well-defined energies in this limit. By 
averaging over time, we force them to slowly sample the whole range of energies between 
— |A| and |A|. Evidently, this would take an infinite amount of time as ujq —>■ 0. 

The Fourier transform G(u) can now be readily determined. In the limit ujq — > 0, 
1/ Juj 2 — |A| 2 with all of the spectral weight 
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A\ 2 -u 2 



(5) 



concentrated between — |A| and |A|. This result can be simply understood. On a time 
scale less than l/u> , each individual fermion has a well-defined frequency u> equal to the 
instanteneous amplitude of the external source. Since the latter fluctuates as |A| cosuot, the 
probability density to find a certain strength uj of this field is given precisely by @. Thus 
the density of fermions whose pole is shifted from to u is (0). Similar broadening in the 
limit of slow fluctuations has been discussed in connection with a pairing || and Peierls @ 
pseudogaps. This simple observation could perhaps vindicate the FLEX approximation in 
light of its recent criticism : the strong broadening of quasiparticle peaks by slow AFM 
fluctuations is a natural phenomenon and not an artifact of the FLEX approximation. 

When ojq ^ 0, the spectral weight of an a± fermion is a superposition of sharp peaks 
A{u) = J2n z n8(u ~ nuJo) with residues z n = J 2 (| A|/o>q). (This result differs markedly from 
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a Poisson distribution found in the case of a quantum boson ||.) The spectral weight at 
uj = is a remnant of the bare fermion. In the fast limit, Zq ~ 1 — \A\ 2 /2u 2 , which agrees 
with the perturbative expansion for z M). In the slow limit, 

^ o;o[l + sin(2|A|/a;o)] 
Z ° ttJA| ' [) 

which certainly does not agree with (£|). Apart from the oscillatory part, z in this limit 
can be obtained as the area under the adiabatic density of states © between —ujq/2 and 
o>o/2, as illustrated in Fig. |^ (note that chaotic oscillations of the peak heights are absent 
close to u — ±|A|). Thus, unlike the perturbation theory, the picture of adiabatic evolution 
provides a reasonable description of the system in the limit of slow fluctuations. 

Fermions with nonzero bare energy: adiabatic limit. When the bare fermion 
modes have unequal energies = =te, the problem does not admit such a simple 
solution for a finite ujq because the Hamiltonian does not commute with itself at different 
times. However, it is still possible to find the r-averaged propagator (0) in the adiabatic 
limit ujq — > 0. To do so, we first determine the propagator for a given instanteneous strength 
of the external source S = |A| cos aj t, 

G{u,5) = 

and then average it over 5 using the probability density (H). As a result, off-diagonal 
(anomalous) components G\2 and G21 vanish, while 




C11H = /_ 



A l dd uj + e / uj + e 



H A I tt^IAI 2 -^ uj 2 -e 2 -5 2 V (u - e)(to 2 - e 2 )' 
with a nonzero spectral weight in the two ranges e 2 < uj 2 < e 2 = e 2 + |A| 2 : 




UJ 



ttv(uj — e)(e 2 — uj 2 ) 



(7) 



For e / 0, there are three peaks in the spectral function, at u — e and ±e. None of these 
is a pole of the propagator, but rather an inverse square-root branch point. In the limit 
e ^> |A|, the cut near e evolves into a simple pole, while the other cut vanishes (Fig. |]). 
At zero temperature, the occupation number of a fermion mode with bare energy e is 



n(e) = [° A u (uj)duj = 1/2 - (e/vre) K'(e/e) 
J — 00 



where K'(&) is a complete elliptic integral of the first kind. There is no discontinuity as 
e crosses zero, although the slope of n(e) diverges logarithmically as e — ■> 0. Assuming a 
constant bare density of fermion states (DOS), it is also possible to obtain the DOS in 
the adiabatic limit by integrating ([!]) over e. The resulting DOS is identical to that of a 
pure ci-wave superconductor, it vanishes linearly with u in the middle of the pseudogap and 
diverges logarithmically at its edges. 
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The simple model considered in this note indicates that the midgap band of Kampf and 
Shrieffer may in fact be inseparable from the two bands above and below the pseudogap. This 
effect can be ascribed to the fluctuations of the gap amplitude in this model. Calculations 
of the FLEX type inevitably contain such fluctuations (see, e.g., |J) and therefore do not 
lead to a sharp energy gap. In principle, this adiabatic broadening of the quasiparticle peaks 
can be lifted if the observation time is shorter than the inverse frequency of the fluctuations 
1/ujq (but still longer than 1/|A|). If fluctuations can also propagate in space, one has to 
limit observations to a spatial area without significant retardation. Given the speed of the 
propagating fluctuations s, the size of such an area should not exceed s/ujq. 

I acknowledge financial support from NEDO (Japan) and the US Department of Energy. 
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FIGURES 



FIG. 1. Feynman diagrams for G(u) through second order in the coupling constant g = |A| 2 /4. 
Broken lines: |A|/2. Solid lines: bare propagators at a frequency shown on the left. 



FIG. 2. Density of states for | A|/ojo = 30. Height of impulses represents spectral weight of 
discrete levels. The broken line is the result of an adiabatic extrapolation. 



FIG. 3. Fermion spectral weight A(u>) in the adiabatic limit for positive bare energies 
= . . . 2|A|. Dotted lines are uj = ±e and u = ±e. 
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FIG. 1. Feynman diagrams for G(uj) through second order in the coupling constant g = |A| 2 /4. 
Broken lines: |A|/2. Solid lines: bare propagators at a frequency shown on the left. 
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FIG. 2. Density of states for |A|/tt>o = 30. Height of impulses represents spectral weight of 
discrete levels. The broken line is the result of an adiabatic extrapolation. 
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